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First results on the measurement of localized fluctuations in the multiplicity of charged particles and photons 
produced in central 158-A GeV/c Pb+Pb collisions are presented. The charged versus neutral correlations in 
common phase space regions of varying azimuthal size are analyzed by two different methods. The measured 
results are compared to those from simulations and to those from different types of mixed events. The compar- 
ison indicates the presence of non-statistical fluctuations in both charged particle and photon multiplicities in 
limited azimuthal regions. However, no correlated charge-neutral fluctuations are observed. 
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The formation of hot and dense matter in high energy 
heavy-ion collisions offers the possibility to create a new 
phase where matter is deconfined and chiral symmetry is re- 
stored. Indications for the formation of such a Quark Gluon 
Plasma (QGP) phase are provided by several results from ex- 
periments at the CERN-SPS [|lj). Event-by-event fluctuations 
in the particle multiplicities and their ratios have recently been 
predicted to provide information about the nature of the QCD 
phase transition Fluctuations may also be caused by 

Bose-Einstein correlations, resonance decays, or more exotic 
phenomena such as pion lasers [Eh. Enhanced fluctuations in 
neutral to charged pions have been predicted as a signature of 
the formation of disoriented chiral condensates (DCC) 
which might be one of the most interesting predicted conse- 
quences of chiral symmetry restoration. 

Theoretical predictions suggest that isospin fluctuations, 
caused by formation of a DCC, would produce clusters of co- 
herent pions in localized phase space regions or domains. The 
probability distribution of the neutral pion fraction in such 
a domain would follow the relation P(f) = 1/2-s/f, where 
/ = N n o /N v . Thus DCC formation in a given domain would 
be associated with large event-by-event fluctuations in the ra- 
tio of neutral to charged pions in that domain. Experimen- 
tally, such fluctuations can be deduced from the measurement 
of fluctuations in the number of photons to charged particles 
in limited r)-<j> regions. The anti-Centauro events, reported by 
the JACEE collaboration with large charged-neutral fluc- 
tuations are possible candidates for DCC events. The studies 
carried out so far in p — p JlO| ] and heavy ion [ pT| , |l2"| | reac- 
tions have searched for fluctuations which extend over a large 
region of phase space. These measurements have provided 
upper limits on the presence of DCC-like fluctuations. 

In this letter we present first results on the search for non- 
statistical event-by-event fluctuations in the relative number of 
charged particles and photons in localized phase space regions 
for central 158 A GeV/c Pb+Pb collisions. The data presented 
here were taken with the 158 -A GeV Pb beam of the CERN 
SPS on a Pb target of 213 /im thickness during a period of 
WA98 operation without magnetic field. The analysis makes 
use of a subset of detectors of the WA98 experiment which 
are used to measure the multiplicities of charged particles and 
photons. Charged particle hits (iV c h) were counted using a cir- 
cular Silicon Pad Multiplicity Detector (SPMD) [[TTJ] located 
32.8 cm downstream from the target. It provided uniform 
pseudorapidity coverage in the region 2.35 < rj < 3.75. The 
detector was 99% efficient for charged particle detection. The 
photon multiplicity was measured using a preshower Photon 
Multiplicity Detector (PMD) |[l3|] placed 21.5 meters down- 
stream of the target and covering the pseudorapidity range 
2.9 < rj < 4.2. It consisted of an array of 53,200 plastic scin- 
tillator pads placed behind 3Xq thick lead converter plates. 
Clusters of hit pads having total energy deposit above a hadron 
rejection threshold are identified as photon— like. The pho- 
ton counting efficiency and the purity of the 7 — like sample 
have been found to be 68% and 65%, respectively, for cen- 
tral events [p4[]. For this analysis the pseudorapidity region 



of common coverage of the SPMD and PMD was selected 
(2.9 < 77 < 3.75). The acceptance in terms of transverse mo- 
mentum (pt) extends down to 30 MeV/c, although no explicit 
PT-selection is applied. Events with pile-up or downstream 
interactions were rejected in the off-line analysis. Strict data 
selection and cleanup cuts have been applied as described in 
Ref. After cuts, a total of 85K central events, corre- 

sponding to the top 5% of the minimum bias cross section as 
determined from the measured total transverse energy, have 
been analyzed. 

The measured results are interpreted by comparison with 
simulated events and with several types of mixed events. Sim- 
ulated events were generated using the VENUS 4.12 jl5] ] 
event generator with default parameters. The output was pro- 
cessed through a WA98 detector simulation package in the 
GEANT 3.21 |16] framework. The centrality selection for 



the simulated data has been made in an identical manner 
to the experimental data by selection on the simulated total 
transverse energy in the WA98 acceptance. The simulated 
VENUS+GEANT events (referred as V+G) were then pro- 
cessed with the same analysis codes as used for the analysis 
of the experimental data. 

The effect of non-statistical DCC-like charged-neutral fluc- 
tuations has been studied within the framework of a simple 
model in which the output of the VENUS event generator has 
been modified. It is expected that DCC domains will occur in 
small regions and will mostly modify the production of low 
momentum pions. The influence of a DCC on the charge- 
neutral pion ratio may then be limited to localized rj-cf> regions 
due to the motion of the DCC domain within the overall col- 
lective motion. To implement the DCC effect, the charges of 
the pions within a localized r\-<\> region predicted by VENUS 
are interchanged pairwise (7r + 7r~ <-> 7r°7r°) according to the 
l/2\/7 probability distribution. The DCC-like fluctuations 
were generated over r\ = 3 — 4 for varying intervals in Acf>. 
Since the probability to produce events with DCC domains is 
unknown, ensembles of events, here referred to as a "nDCC 
events", were produced which consisted of a mixture of nor- 
mal events with varying fractions of pure DCC-like events. 
The nDCC events were then tracked through GEANT. 

In the search for evidence of non-statistical charged-neutral 
fluctuations, two different analysis techniques have been ap- 
plied. The first method employed in the present analysis is the 
technique of discrete wavelet transformations (DWT). DWT 
methods are now widely used in many applications, such as 
data compression and image processing, and have been shown 
to provide a powerful means to search for localized domains 
of DCC Jl7| , |l8| ]. While there are several families of wavelet 
bases distinguished by the number of coefficients and the 
level of iteration, we have used the frequently employed D— 4 
wavelet basis [20], which are orthogonal, continuously dif- 
ferentiable and localized in space. The analysis has been per- 
formed with the sample function chosen to be the photon frac- 
tion, given by /'((/>) = N- ( -rM4>)/(N-y- mc ((t>) + N ch (4>)) 
as a function of the azimuthal angle <f>, with highest resolu- 
tion scale j m ax = 5. The input to the DWT analysis is the 
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FIG. 1. The rms deviations of the FFC distributions at j = 1 
for simulated nDCC events with extent Atfincc = 90° and for vari- 
ous mixed events constructed from those events, as a function of the 
fraction of DCC-like events present in the nDCC sample. 



spectrum of the sample function at the smallest bin size cor- 
responding to the highest resolution scale, j m ax, where the 
number of bins is 2 Jmax . The sample function is then ana- 
lyzed at different scales j by being re-binned into 2 3 bins. 
The DWT analysis yields a set of wavelets or father function 
coefficients (FFC) at each scale from j = l,...,(j maa; — 1). The 
coefficients obtained at a given scale, j, are derived from the 
distribution of the sample function at one higher scale, j + 1. 
The FFCs quantify the deviation of the bin-to-bin fluctuations 
in the sample function at that higher scale relative to the av- 
erage behavior. The presence of localized non-statistical fluc- 
tuations will increase the rms deviation of the distribution of 
FFCs and may result in non-Gaussian tails jl7| , ^8[ ]. 

The DWT technique as used in this analysis is demonstrated 
in Fig. [I] where it has been applied to the simulated nDCC 
events. The rms deviation of the FFC distribution is shown 
as a function of the fraction of DCC-like events in the nDCC 
sample. The rms deviation is observed to increase strongly 
with increasing DCC-like fraction. Due to the inherent un- 
certainties in the description of "normal" physics and detector 
response in the V+G simulations, the observation of an exper- 
imental result with rms which differs from the case with zero 
DCC fraction cannot be taken alone as evidence of DCC ob- 
servation. For this reason four different types of mixed events 
have been created from the real or simulated events in order 
to search for non-statistical fluctuations by removing various 
correlations in a controlled manner while preserving the char- 
acteristics of the measured distributions as accurately as pos- 
sible. The first type of mixed events (Ml), are generated by 
mixing hits in both the PMD and SPMD separately, with no 
two hits taken from the same event. Hits within a detector in 
the mixed events are not allowed to lie within the two track 
resolution of that detector. The second kind of mixed events 
(M2) are generated by mixing the unaltered PMD hits of one 
event with the unaltered SPMD hits of a different event. In- 
termediate between the Ml and M2 kinds of mixed events is 
the case where the hits within the PMD are unaltered while 



the SPMD hits are mixed (M3-7), or the SPMD hits are un- 
altered while the PMD hits are mixed (M3-ch). In each type 
of mixed event the global (bin 1) Ar 7 _iike— A^h correlation is 
maintained as in the real event. 

The rms deviations of the FFCs for the different kinds of 
mixed events produced from the nDCC events are also shown 
in Fig. [j]. In the case of vanishing DCC-like fluctuations, the 
rms values of the various types of mixed events are very close 
to each other. The V+G rms values are lower than those of the 
mixed events due to the presence of additional correlations 
between and -/V 7 _ijke mostly as a result of the charged 
particle contamination in the A^-ijke sample. The rms devi- 
ations for the Ml events are found to be almost independent 
of probability of DCC-like events, while the rms deviations of 
the M2 events increase similarly, but more weakly, than those 
of the nDCC events. The rms deviations for the M3 sets of 
events are found to lie between M2 and Ml. Thus, the se- 
quence of the mixed events relative to the simulated events 
(or data) gives a model independent indication of the presence 
and source of non-statistical fluctuations. The simple DCC 
model used here results in an anti-correlation between iV 7 _iiko 
and A^ch due to the "isospin-flip" procedure used to implement 
the DCC effect. It also results in non-statistical fluctuations in 
both iV 7 _iikG and N c ^. Thus the M2 events remove only the 
Af-y-iike-A'ch anti-correlation while the Ml events are seen to 
remove all non-statistical fluctuations and correlations. The 
M3 mixed events give intermediate results because they con- 
tain only the iV 7 _iik (M3-7) or N c ^ (M3-ch) non-statistical 
fluctuations. 

The FFC distributions extracted from the measured f'(4>) 
ratio are shown in the bottom panel of Fig. || for the experi- 
mental data, for Ml events (from data), and for V+G events. 
The results are shown for scales j —I and 2, which carry in- 
formation about fluctuations at 90° and 45° in (j). The FFC 
distributions of the experimental distributions are seen to be 
broader than the V+G and Ml results. This suggests the pres- 
ence of non-statistical fluctuations. 

A more conventional method similar to that described in 
Ref. [nljl has also been used to search for non-statistical fluc- 



tuations. The correlation between N. 



7— like 



and A^ch has been 



studied in varying <j) intervals, by dividing the entire 0-space 
into 2, 4, 8, and 16 bins. The correlation plot of A^h ver- 
sus A r 7 _iik c [ 19] is obtained for each (f> segmentation, starting 
with the case of 1 bin which corresponds to the full azimuth. 
A common correlation axis (Z) has been obtained for the full 
distribution by fitting the mean A^ 7 _iik and A'ch values with a 
second order polynomial. The distance of separation (Dz) be- 
tween the data points and the correlation axis has been calcu- 
lated with the convention that Dz is positive for points below 
the Z-axis (increasing N^-\^ c ). The distribution of Dz rep- 
resents the fluctuations of A^-iikc relative to A^h compared 
to the common correlation axis. In order to compare fluctu- 
ations at different bin sizes having different multiplicities we 
use a scaled variable, Sz = Dz / s(Dz), where s(Dz) is the 
rms deviation of the Dz distributions for V+G events. The 
presence of events with localized non-statistical fluctuations 



3 



10 5 










4 


: 4 bins jr 




: 8 bins jf 




10 






/ 




10 










10 2 










10 










1 


B.,.Ih..,, 






,,;,L, 



-10 





10 5 






ent 


10 4 


> 






10 3 


O 


10 2 


u 






1 


10 






z 


1 



-5 5 



-5 5 10 



r 4 bins 




r 8bimi^^ 








r f 

; . . . . ii '. . . . i . . 


L. 





-0.5 0.5 ,-,,-,„ -0.5 0.5 1 



FIG. 2. The Sz and FFC distributions for 4 and 8 divisions in 
cj>. The experimental data, Ml and V+G events are shown by solid 
circles, solid histograms and dashed histograms, respectively. The 
statistics for data and mixed events are same, whereas the distribution 
for the V+G events is normalized to the number of data events. 
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would be expected to result in a broader distribution of Sz 
compared to those for normal events. The Sz distributions 
calculated at 4 and 8 bins in <f> angle are shown in the top panel 
of Fig. U for data, Ml, and V+G events. The experimental 
distributions are broader than the simulation and Ml results, 
again indicating the presence of additional fluctuations. 

The rms deviations of the Sz and FFC distributions as a 
function of the number of bins in azimuth is shown for exper- 
imental data, mixed events, and V+G in Fig. ^. The statistical 
errors on the values are small and lie within the size of the 
symbols. The error bars include both statistical and systematic 
errors. The systematic errors have been estimated by investi- 
gation of effects such as the uncertainties in the detection ef- 
ficiencies, gain fluctuations, backgrounds, binning variations, 
and fitting procedures. Since the mixed events are constructed 
to maintain the iV 7 _uke— correlations for the full azimuth 
(bin 1), the rms deviations of data and mixed events for this 
bin are identical. The difference of the Sz rms deviations be- 
tween data and V+G for this bin is the same as reported in 
an earlier WA98 publication JTl|]. The comparison of V+G 
and the Ml mixed events demonstrates the utility of the DWT 
method to normalize out the average behavior when the bin- 
to-bin fluctuation information is extracted. 

As noted in the discussion of Fig. Ill even in the absence of 
DCCs there exist uninteresting correlations between iV 7 _iik c 
and iV c h which are removed by the event mixing procedure 
and thereby result in a difference between the real and mixed 
events. The various mixed event rms values of Fig. || have 
therefore been rescaled by the percentage difference between 
the rms deviations of the V+G distributions and those of the 
corresponding V+G mixed events in order to better illustrate 



effects in the data beyond those present in V+G. Taking the 
larger of the asymmetric systematic error at each point to be 
one sigma, we find that at 2, 4, and 8 bins the values of the 
Sz rms deviations of the data are 3.2a, 3.4er, and 3.1er larger 
than those of Ml events, respectively. Similarly, the FFC rms 
deviations at 4 and 8 bins for data are 4.3cx and 3.3cr larger 
than those of the Ml events. At 16 and 32 bins, the result 
for mixed events and data agree within the quoted errors. The 
corrected rms deviations of the M2 events agree with those of 
the experimental data within error for all bins. The M3 type 
mixed events are found to be similar to each other within the 
quoted errors and lie between Ml and M2. 

The observation that the rms deviations of the Sz and FFC 
distributions for experimental data are larger than those of the 
Ml events indicates the presence of localized non-statistical 
fluctuations. However, the comparision of the rms deviations 
for data with those of M2 events implies the absence of event- 
by-event correlated fluctuations in N^-jxke versus N^. The 
M3-type mixed events indicate the presence of localized inde- 
pendent fluctuations in iV 7 _iikc and iV c h of similar magnitude. 

If the amount of DCC-like fluctuations in the experimental 
data were large, then the rms deviations shown in Fig. || for 
data would have been larger compared to those of M2 events. 
Since this is not the case, we compare the measured results 
with those obtained from the simulation as shown in Fig. [l] 
to extract upper limits on the probability of DCC-like fluctu- 
ations at the 90% confidence level. Within the context of this 
simple DCC model, upper limits on the presence of localized 
non-statistical DCC-like fluctuations of 10~ 2 for Acf> between 
45-90° and 3 x 10~ 3 for A<f> between 90-135° are extracted. 

In summary, a detailed event-by-event analysis of the fluc- 
tuations in the rj — <f> phase space distributions of charged 
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particles and photons has been performed for central Pb+Pb 
collisions at 158-^4 GeV using two complementary analysis 
methods. The first analysis employed the discrete wavelet 
transformation technique to investigate the relative magni- 
tude of the A^ 7 _iiko versus iV c h fluctuations in adjacent phase 
space regions. The second method studied the magnitude 
of the TV-y-Hke versus iV c h multiplicity fluctuations in de- 
creasing phase space regions. The results were compared to 
pure VENUS+GEANT simulations and to various types of 
mixed events to isolate the source of non-statistical fluctua- 
tions. Both analysis methods indicated non-statistical fluctua- 
tions beyond simulation and beyond pure mixed events at the 
3-4ct level for (j> intervals of greater than 45°. This is found to 
be due to uncorrected non-statistical fluctuations in JV 7 _iik e 
and iV c h. No significant correlated fluctuations in iV 7 _uko ver- 
sus iV c h were observed. The results allow to set an upper limit 
on the frequency of production of DCCs of limited domain 
size, as demonstrated with a simple model of DCC-like fluc- 
tuations. 
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